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Summary:  An  experimental  investigation  of  the  evolution  of  the  profiles  of 
a  set  of  spilling  and  a  set  of  plunging  breaking  waves  is  presented.  Each 
set  of  waves  consists  of  similar  waves  scaled  to  three  characteristic  wave 
periods.  Measurements  of  various  wave  heights,  periods,  and  phase  speeds 
are  presented. 

1.  Introduction \ 

Breaking  waves  are  important  in  many  oceanographic  phenomena  and  in  the 
design  of  marine  structures  and  vehicles.  In  oceanography,  wave  breaking 
provides  a  mechanism  In  which  wind  generated  waves  transfer  their  energy 
to  turbulence.  This  turbulence  enhances  the  transfer  of  heat  and  gases 
across  the  air-sea  Interface.  The  forces  produced  by  breaking  waves  on 
ships  and  marine  structures  can  also  be  of  great  Importance.  At  the 
instant  that  the  wave  breaks,  the  fluid  velocities  near  the  crest  are 
about  equal  to  the  wave  phase  velocity.  The  momentum  of  this  mass  of 
water  can  produce  Impact  pressures  that  are  capable  of  capsizing  boats  and 
damaging  marine  structures. 

In  order  to  study  any  of  the  above  phenomena  in  the  laboratory,  one  must 
first  produce  a  set  of  waves  that  cover  an  appropriate  range  of  sizes  and 
types  to  drive  the  phenomena  under  study.  This  is  not  a  simple  task 
because  breaking  waves  are  not  well  understood.  Fundamental  to  this 

problem  is  the  Identification  of  a  simple  set  of  Independent  variables 
that  quantitatively  describe  a  breaking  wave.  In  the  case  of  steady, 
periodic,  nonbreaking  waves,  the  Independent  variables  have  been  known 
since  the  original  work  of  Stokes  (1847).  Here  the  wave  is  completely 
described  by  two  Independent  variables,  either  its  period  and  amplitude, 
or  its  wavelength  and  amplitude.  In  the  breaking-wave  case,  the  wave 
profile  evolves  as  it  propagates  changing  in  period  and  amplitude  until 
the  crest  either  spills  or  plunges  forward.  This  evolution  can  be  driven 
by  varying  water  depth,  Interaction  with  a  mean  current,  or,  as  in  the 
present  experiments,  interaction  with  other  components  in  a  wave  packet. 
The  variables  that  completely  describe  these  waves  are  not  known.  It 
seems  reasonable  that  one  of  these  would  be  the  apparent  period  of  the 
wave  at  the  time  it  begins  to  break.  As  in  the  steady  wave  case,  the 
period  should  be  related  to  the  wavelength  (but  probably  in  some  complex 
way),  thus  determining  the  general  length  scale  of  the  wave.  For  a  wave 
with  a  given  period  at  the  point  of  incipient  breaking,  it  is  thought  that 
one  or  more  variables  will  determine  the  breaker  type  which  can  range  from 
a  gentle  splller  to  a  plunging  breaker.  These  variables  might  include 
wave  slope  and  asymmetry  parameters  at  the  time  of  breaking,  as  well  as 
parameters  Chat  depend  on  the  time  history  of  the  wave  profile  (for 
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example,  the  rate  of  lacrease  of  wave  height  or  steepness). 

There  are  two  goals  in  the  present  experimental  research  program.  First, 
we  wish  to  develop  a  technique  for  generating  a  given  breaker  type  at 
various  wave  periods  with  a  single  wavemaker.  These  waves  could  then  be 
used  In  subsequent  studies  to  test,  for  Instance,  the  capsize  resistance 
of  boat  models  In  the  presence  of  several  plunging  and  spilling  breakers 
as  a  function  of  wave  period.  Our  second  goal  Is  to  search  for  a  set  of 
Independent  variables  that  describe  these  waves  and,  In  particular,  to 
find  quantitative  measures  that  Indicate  the  occurrence  of  breaking  and 
the  breaker  type.  Several  papers  In  which  the  profiles  of  breaking  waves 
have  been  studied  quantitatively  have  appeared  In  the  recent  literature 
(Bonmarln  and  Ramamonjlarlsoa  ,  1984,  KJeldsen  and  Myrhaug,  1979,  and 

Melville  and  Rapp,  1985).  In  the  present  study  we  expand  upon  this 
previous  work  by  examining  one  spilling  and  one  plunging  breaking  wave, 
each  of  which  Is  scaled  to  three  absolute  wave  periods. 

2.  Experimental  Details. 

2a.  Facilities. 

The  experimental  research  program  was  conducted  In  a  wave  tank  at  the 
United  States  Naval  Academy.  The  tank  Is  36.6  m  In  length,  2.44  m  wide  and 
1.5  m  deep.  It  Is  equipped  with  an  MTS  dual  flap  wavemaker  which  Is  2.44 
m  wide  and  1.83  m  high,  extending  the  full  tank  depth.  A  hydraulic 
actuator  Is  used  to  drive  the  lower  wave  board  with  respect  to  the  tank 
foundation;  a  second  actuator  Is  used  to  drive  the  upper  wave  board  with 
respect  to  the  lower  wave  board. 

2b.  Wave  Height  Measurments. 

Wave  profile  measurements  were  obtained  with  four  linear  resistance  wave 
height  gauges.  The  gauges  were  mounted  20  to  30  cm  apart  along  the  tank 
centerline  on  a  movable  platform.  For  each  wave,  measurements  were  taken 
with  the  platform  at  5  separate  positions,  resulting  In  wave  records  at  14 
locations.  Positions  were  chosen  ranging  from  a  small  distance  past  the 
break  point  to  about  one  wavelength  toward  the  wavemaker  from  the  break 
point.  The  heights  of  the  wave  crest  and  troughs,  and  the  various  wave 
periods  were  measured  from  each  single  wave  record.  Us'lng  the  probe 
pairs,  the  speeds  of  the  crest,  troughs,  and  zero  crossings  were 
computed.  Measurements  were  repeated  three  times;  reported  values  are  the 
average  of  these  results. 

3.  Wavemaker  Driving  Signal 

The  computer  program  that  generated  the  driving  signal  for  the  wavemaker 
was  a  modification  of  the  program  used  by  Salslch  et  al  (1983).  The 
purpose  of  the  modifications  was  to  facilitate  the  scaling  of  a  given 
driving  signal  to  various  frequencies.  In  Its  modified  form,  the  program 
generates  a  signal  that  can  be  divided  Into  two  basic  segments.  In  the 
first  segment,  both  the  frequency  and  amplitude  of  the  signal  were  held 
constant.  During  this  period,  the  wavemaker  hydraulic  system  was  brought 
up  to  full  stroke  while  producing  a  few  very  small  amplitude  waves .  In 
the  second  segment,  the  frequency  was  decreased  linearly  with  time.  The 
sweep  In  the  frequency  during  this  segment  causes  the  wave  pattern  to 
converge^  since  the  lower  frequency  waves  overtake  the  slower  waves 
generated  ahead  of  them.  By  linear  theory  (Longue t-Hlgglns  1974)  the 


distance  from  the  wavemaker  to  the  point  where  the  waves  converge  is  given 
by 
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Att  df/dt 

where  f  is  the  frequency  in  hertz,  t  is  time  and  g  is  the  acceleration  of 
gravity*  Wave  breaking  occurs  at  scnne  distance  from  the  wavemaker  which 
is  less  than  X^,.  The  amplitude  of  the  signal  is  also  varied  as  the 
frequency  is  varied.  As  the  frequency  decreases  from  f ,  its  initial 
value,  to  an  intermediate  frequency,  the  amplitude  is  increased.  Then, 
as  the  frequency  continues  to  decrease  to  £2 ,  the  amplitude  decreases  to 
its  original  value.  In  order  to  describe  this  amplitude  envelope,  we 
first  define  a  frequency,  f^,  during  the  time  when  the  frequency  is 
between  fi  and  f„: 

^  (f,  -  f 
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Note  that 
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where  f^^  ■  0.5(fj  +  f2). 

The  amplitude  of  the  wave  envelope  during  this  phase  of  the  signal  is  then 
given  by 
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where  E  and  P  are  positive  constants, 
another  frequency,  f^,  is  defined: 

^^2  -  'av> 


=  ^2 


Note  that 
f d  ■ 

and 


fav  at  f 


Between  frequencies  f  and  f2 
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fjj  “  f2  at  f  “  f2 . 

The  amplitude  of  the  wave  envelope  in  this  region  is  given  by 
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The  amplitude  functions  have  the  properties:  A^Cfi)  ■  Ajj(f2)»  and  “ 

Ad  ^  “  1.0. 

In  the.study  described  below,  we  have  chosen  a  set  of  vedues  for  E,  P, 
fp/fl,  and  f2/fl  and  scaled  the  wave  generation  signal  to  several  values 
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of  absolute  frequency!  say  f  •  Given  the  above  parameters}  we  then  adjust 
the  theoretical  break  pointi^X^^,  in  each  case  to  create  thq  same  number  of 
waves  in  the  generation  signal.  Let  the  variables  describing  two  scaled 
wave  signals  be  distinguished  by  *  and  The  first  step  in  making  the 
signals  similar  is  to  scale  the  three  frequencies  (f^  f^,  and  by  the 
same  ratio: 


f  "/f 
1  '  1 


7f  '  -  fo"/f  '  -a. 

p  2  2 


Since  the  frequency  is  the  number  of  waves  passing  a  given  point  per  unit 
time,  the  nvunber  of  waves  generated,  n,  is  given  by  the  integral ‘of  the 
frequency  over  time.  Thus, 


where  df/dt  is  given  by  Equation  (1)  and  is  obtained  from 

^2  ■  df 

T~o  -  dF 

After  performing  the  integration  and  some  algebraic  manipulation.  Equation 
(7)  becomes 

27r  X, 

n - -  (f,^-  f/)  (9) 

g  1  2 

Since  the  number  of  waves  in  the  two  wave  trains  is  to  be  the  same  we  have 
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Or  using  Equation  (6) 


Because  the  transfer  function  between  the  wavemaker  and  the  waves 
generated  is  a  nonlinear  function  of  frequency,  one  would  not  expect  the 
wave trains  generated  by  these  scaled  driving  signals  to  have  exactly  the 
same  scaled  profiles.  However,  due  to  the  preliminary  nature  of  this 
study,  no  attempt  to  account  for  these  nonllnearltles  was  made.  The 
results  presented  below  indicate  that  within  the  range  of  frequencies 
used,  the  scaling  procedure  was  adequate. 

4. Results  and  Discussion 


The  waves  used  in  this  study  were  produced  with  one  tyne  of  wavemaker 
driving  signal.  The  parameters  for  this  signal  were  fp/f^^O.TO, 
f2/fi"0.60,  E  “  5.0,  and  P  “  4.0.  A  plot  of  the  signal  appears  in  Figure 
1.  This  signal  was  scaled  to  three  different  frequencies  (fp  “0.70,  0.85, 
and  1.0).  For  each  frequency,  two  amplitudes  of  the  signal  were  chosen  so 


that  the  wave  broke  at  x/X>,“  0.67  and  0.80,  where 


is  the  distance  from 


the  wavemaker;  the  breaking  events  at  x/Xb"  0*80  were  spilling  breakers, 
while  the  breaking  events  at  x/Xb“  0.67  were  plunging  breakers.  These 


results  are  consistent  with  Salslch  et  al  (1983).  When  the  signals  were 
first  Input  to  the  wavemaker,  two  breaking  events  sometimes  occurred  near 
the  desired  point  of  Incipient  breaking  (break  point).  This  problem  was 
remedied  by  varying  the  phase  of  the  drive  signal  within  the  envelope 
described  by  equations  (3),  and  (5)  until  a  single  breaking  event 
occurred . 

A  set  of  wave  profiles  taken  at  four  positions  near  the  break  point  for 
one  of  the  waves  appears  In  Figure  2.  Wave  profiles  such  as  these, 
measured  within  one  wavelength  of  the  break  point,  were  used  to  study  the 
behavior  of  the  wave  as  It  approached  breaking. 

The  most  fundamental  parameter  of  Interest  Is  the  wave  period  measured  by 
a  single  probe.  A  sketch  showing  the  definitions  of  the  various  features 
of  the  breaking  wave  profile  Is  given  In  Figure  3.  The  wave  period,  T^j.ls 
defined  as  the  time  between  the  passage  of  the  downward  zero-crossing  of 
the  trough  ahead  of  the  breaker  and  the  downward  zero-crossing  of  the 
breaking  crest  (T^j*  T^  -  T^).  A  second  wave  period,  called  the  crest 
period,  T^ ,  la  defined  as  the  time  between  the  passage  of  the  upward 
zero-crossing  ahead  of  the  breaking  crest  and  the  downward  zero-crossing 
after  the  crest  (T  ■  Tj  -  T^).  Figure  4a  contains  a  plot  of  these  wave 
periods  versus  distance,  z,  from  the  nominal  break  point  for  the  three 
spilling  breakers,  and  Figure  4b  contains  the  data  for  the  three  plunging 
breakers.  The  horizontal  distance,  x.  Is  scaled  by  the  wavelength,  L^,  of 
a  linear  wave  whose  period,  T^,  equals  the  value  of  T^  at  the  break  point 
In  each  case.  The  break  point,  which  Is  near  x  *  0,  was  chosen  based  on 
the  plots  of  crest  height  versus  distance  from  the  break  point  (see 
below).  The  periods  are  nondlmenslonallzed  by  the  appropriate  T^.  Notice 
that  In  all  six  cases  both  periods  decrease  slowly  as  the  break  point  Is 
approached.  In  the  vicinity  of  the  break  point,  T^  reaches  a  minimum  and 
then  Increases,  while  T^  becomes  nearly  constant.  The  significance  of 
this  minimum  Is  not  known.  The  collapse  of  the  data  with  the  above 
scaling  Is  excellent  for  the  spilling  breakers  and  good  for  the  plunging 
breakers. 

The  height  of  the  crests  of  the  breakers,  »  and  the  depths  of  the 
preceding  and  the  following  troughs,  n^-p  and  n^f,  are  plotted  versus 
horizontal  distance  from  the  break  point  In  Figure  5a  for  the  spilling 
breakers  and  In  Figure  5b  for  the  plunging  breakers.  All  variables  are 
nondlmenslonallzed  by  the  appropriate  value  of  for  each  wave.  The 
height  of  the  crest  Increases  as  the  break  point  Is  approached.  In  all 
but  one  case,  the  plunging  breaker  with  f  ■  1.0,  the  crest  height  reaches 
a  maximum  at  the  point  where  the  wave  vlsuUly  appears  to  break.  This  Is 
In  agreement  with  the  experimental  results  of  Bonmarln  and  Ramamonjlarlsoa 
(1984)  and  the  numerical  results  of  Longuet-Hlgglns  and  Cokelet  (1976). 
The  point  where  the  height  of  the  crest  reached  a  maximum  was  chosen  as 
the  point  of  Incipient  breaking  (break  point)  In  those  five  cases.  In  the 
sixth  case,  the  plunging  breaker  with  fp  “  1.0,  the  data  point  with  the 
smallest  x  was  taken  as  the  break  point.  Thus,  the  break  point  was 
between  0  and  0.1  x/Lq  In  all  cases.  With  the  exception  of  the  plunging 
breaker  with  fp  ■  1.0,  the  data  collapses  quite  well.  The  dimensionless 
wave  height,  /L  ,  at  the  break  point  for  the  spilling  breakers  was  about 
0.08,  while  for  the  Lw^  scaled  plunging  breakers  It  was  0.09.  In  much  of 
the  data,  the  results  for  the  plunging  breaker  with  f  ■  1.0  differ  from 
the  other  two  plungers.  The  reason  for  this  discrepancy  Is  assumed  to  be 
experimental  error,  however,  there  was  Insufficient  time  prior  to  the 
present 'publication  date  to  explore  this  problem  further. 


As  the  break  point  is  approached,  the  depth  of  the  trough  preceding  the 
breaking  crest  decreases,  while  the  opposite  occurs  for  the  following 
trough.  Bonmarln  and  Ramamon j iarlsoa  (1984)  report  a  qualitatively 
similar  behavior  for  the  preceding  trough  in  their  experiments.  The 
curves  for  the  preceding  and  following  troughs  cross  at  x/Lq  “0.1  for  the 
spilling  breakers  and  x/Lq  "  0.35  for  the  plunging  breakers.  The  total 
wave  height  based  on  the  preceding  trough  Increases  slightly  as  the  break 
point  Is  approached,  reaching  0.11  and  0.12  Lq  for  the  spilling  and 

plunging  breakers,  respectively. 

Three  dimensionless  parameters  that  quantify  the  asymmetry  of  the  wave 
form  (sqi  ,  pjj,  and  see  lAHR  1986)  are  plotted  versus  x/L^  In  Figure  6a 
for  the  spilling  breakers  and  In  Figure  6b  for  the  plunging  breakers  (see 
Figure  3  for  definitions).  These  parameters  were  first  used  by  Kjeldsen 
and  Myrhaug  (1979)  who  called  them  e,  p,  and  x»  respectively.  The  crest 
front  steepness,  s^^i,  can  be  considered  as  the  ratio  of  the  average  rise 
velocity  of  the  water  surface  (from  the  time  of  passage  of  the  upward 
zero-crossing  ahead  of  the  crest  to  the  time  of  passage  of  the  crest), 
divided  by  the  celerity  of  a  linear  wave  based  on  the  period,  T^j .  The 
parameter  S(;>  Increases  slowly  up  to  the  break  point  and  then  decrease  for 
all  but  the  plunger  with  fp  ■  1.0.  The  average  value  at  the  break  point  is 
about  0.50  for  the  spilling  breakers  and  0.70  for  the  plunging  breakers. 
Kjeldsen  and  Myrhaug  (1979)  found  Sqi  In  the  range  0.32  to  0.78  In  their 
laboratory  data,  while  Bonmarln  and  Ramamonjlarsoa  (1984)  found  a  maximum 
of  0.55  In  their  Investigation  of  a  single  wave. 

The  asymmetry  parameter  pu  Is  the  ratio  of  the  crest  height,  >  to  the 
vertical  distance  between  tne  preceding  trough  and  the  crest,  ri^.  “  np^- • 
It  Is  a  measure  of  the  asymmetry  of  the  profile  about  the  horlzoiftal 
axis.  This  parameter  continually  Increases  for  all  waves  and  has  an 
average  value  of  about  of  0.75  for  both  the  splllers  and  the  plungers  at 
the  break  point.  Thus,  for  the  waves  studied,  It  seems  that  a  value  of  Wji 
near  0.75  may  indicate  the  occurence  of  a  breaking  event,  while  the  value 
of  Sj^,  at  that  Instant  may  Indicate  the  breaker  type. 

The  last  asymmetry  parameter,  p  ,  Is  a  measure  of  the  asymmetry  of  the 
profile  of  the  wave  crest  aboul^  a  vertical  axis.  It  Is  defined  as  the 
ratio  of  the  time  delay  between  the  passage  of  the  crest  and  the  following 
downward  zero-crossing  to  the  time  delay  between  the  preceding  upward 
zero-crossing  and  the  crest.  This  parameter  Increases  at  first  ard 
reaches  a  maximum  a  little  ahead  of  the  break  point  in  all  cases  except 
the  plunging  breaker  with  f  ■  1.0.  The  average  peak  value  Is  about  1.9  for 
the  plunging  breakers  ana  1.8  for  the  spilling  breakers.  At  the  break 
point,  has  decreased  to  1.7  and  1.3  for  the  plunging  and  spilling 

breakers,  respectively. 

Plots  of  the  horizontal  speed  of  the  crest,  ,  divided  by  the  speed  of 
the  upward  zero-crossing  preceding  the  crest,  C3 ,  and  divided  by  the 
speed  of  the  downward  zero-crossing  following  the  crest,  C5 ,  versus 
horizontal  position  are  given  In  Figures  7a  and  7b  for  the  spilling  and 
plunging  breakers,  respectively.  The  speeds  were  obtained  by  measuring 
the  time  delay  In  the  passage  of  the  feature  In  question  at  two  wave 
height  probes  with  a  known  distance  between  them.  These  phase  speed 
ratios  Indicate  the  average  rate  of  steepening  of  the  forward  and  rear 
face  of  the  wave  crest.  When  c^/ct  >  1,  the  forward  face  of  the  wave  Is 
steepening,  while  the  rear  face  of  the  wave  steepens  when  c^/c^  <  1.  Note 
from  Figure  2  that  the  average  steepness  or  slope  Is  considerably  less 


than  the  Instantaneous  slope  near  the  crest  just  prior  to  breaking.  The 
curves  in  Figure  7  are  jagged  due  to  the  magnification  of  small  errors 
when  subtracting  the  time  of  passage  of  a  given  feature  at  neighboring 
probes;  however,  some  general  trends  can  be  seen  in  the  data.  For  both 
the  spilling  and  plunging  breakers,  the  forward  and  rear  faces  of  the  wave 
are  steepening  slowly  for  large  x  1.0,  c^/ccT  1.0).  As  x 
decreases,  decreases  steadily;  concurrently,  ^c'^3  first 
Increases  and  then  decreases.  At  the  break  point,  and  0^/03  are 
nearly  equal  with  a  value  of  0.8  and  0.9  for  the  spilling  and  plunging 
breakers,  respectively. 

5.  Conclusion 

t 

The  collapse  of  the  data  for  the  three  spilling  and  the  three  plunging 
breakers  indicates  that  the  proposed  method  is  capable  of  generating  a 
given  breaker  type  at  a  range  of  wave  periods.  In  future  work,  attempts 
will  be  made  to  improve  the  accuracy  of  the  method  by  accounting  for  the 
nonlinearity  of  the  transfer  function  of  the  wavemaker.  The  analysis  of 
various  wave  profile  parameters  measured  with  single  wave  height  probes 
indicates  that  the  horizontal  asymmetry  factor,  may  be  the  best 
indicator  of  the  occurrence  of  a  breaking  wave.  It  Increased  steadily  to 
about  0.75  at  the  point  of  incipient  breaking  for  the  six  waves 
Investigated  herein.  During  the  beginning  of  the  breaking  process  it 
continued  to  Increase.  The  crest  front  steepness,  s^, ,  also  Increased  as 
the  wave  approached  the  point  of  incipient  breaking.  It  then  decreased 
during  the  beginning  of  the  breaking  process.  Its  value  at  the  break 
point  was  0.50  for  the  spilling  breakers  and  0.70  for  the  plunging 
breakers.  Thus,  the  value  of  s^i ,  at  the  Instant  when  ■  0.75  may  be  a 
good  indicator  of  breaker  type.  The  use  of  two-probe  arrays  for  detecting 
breaking  events  from  wave  phase  speeds  will  require  Improvements  in  the 
accuracy  of  the  measurements. 
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